Pedestrian Approach to the Two-Channel Kondo Model 



P. ColemanS L. loffe^ and A. M. Tsvelik^ 
^ Serin Laboratory, Rutgers University, P.O. Box 849, Piscataway, New Jersey 08855-0849 
^Department of Physics, University of Oxford, 1 Kehle Road, Oxford 0X1, UK 

We reformulate the two-channel Kondo model to explicitly remove the unscattered charge degrees 
of freedom. This procedure permits us to move the non- Fermi liquid fixed point to infinite coupling 
where we can apply a perturbative strong-coupling expansion. The fixed point Hamiltonian involves 
a three-body Majorana zero mode whose scattering effects give rise to marginal self-energies. The 
compactified model is the N = Z member of a family of "0(iV)" Kondo models that can be solved 
by semiclassical methods in the large N limit. For odd N , fermionic "Kink" fluctuations about the 
N = oo mean-field theory generate a fermionic iV-body bound-state which asymptotically decouples 
at low energies. For N = Z, our semi-classical methods fully recover the non-Fermi liquid physics 
of the original two channel model. Using the same methods, we find that the corresponding 0(3) 
Kondo lattice model develops a spin-gap and a gapless band of coherently propagating three-body 
bound-states. Its strong-coupling limit offers a rather interesting realization of marginal Fermi 
liquid. 
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I. INTRODUCTION 

The possible origins of non-Fermi liquid behavior in 
highly correlated metals is a topic of lively debate. 
Non-Fermi liquid behavior has been observed in a va- 
riety of metallic systems, including low dimensional 
conductors,^ cuprate superconductors,®'^ systems at a 
quantum critical point® and in certain heavy fermion 
metals. ^'"""^ In this context, the two channel Kondo model 
is of particular interest, for it provides a simple exam- 
ple of a class of non- Fermi liquid (NFL) behavior that is 
driven by local physics. ■'■^^■'■^ Variants of this model are of 
potential relevance to the physics of dilute heavy fermion 
systems^' ■'■^'■'■^ and two level tunneling systems. -"-^ 

Most pedestrian techniques in condensed matter the- 
ory involve perturbative expansions about weak or strong 
coupling fixed points. Insight into the underlying physics 
is often gained when these simple approaches can be ap- 
plied. Unfortunately, non-Fermi liquid behavior in the 
two-channel Kondo model is associated with an interme- 
diate coupling fixed point, placing it beyond the reach 
of elementary methods. In this paper show how this diffi- 
culty can be overcome by taking advantage of spin-charge 
decoupling. A key result is the derivation of an effective 
Hamiltonian for the low-energy spin physics of the two- 
channel Kondo model: 

H* = J2 ek*k • *k + «**^^H0)*(^)(0)*(^)(0). (1) 

k 

Here * = *(2), *(3)) 

is a three-component Majo- 
rana fermion (\I>('') = vj/C")*) with dispersion eg and $ = 
is a Majorana fermion localized at the origin, whose 
singular scattering effects generate a marginal Fermi liq- 
uid. The spin density is given by S(x) = — i\f (x) x \f (x). 



We shall derive this Hamiltonian in two ways and show 
how it may be used to simply obtain the thermodynamics 
and spin-physics of the two-channel Kondo model. 

Spin-charge decoupling, the separation of spin and 
charge into independent degrees of freedom at low en- 
ergies, is a universal feature of impurity Kondo models. 
For multi-channel Kondo models, the use of a lattice cut- 
off interferes with this separation. By compactifying the 
two-channel model into a Hamiltonian that exclusively 
describes the decoupled spin degrees of freedom, we show 
how the over-screened fixed point-*--*- is eliminated, relocat- 
ing the NFL fixed point to infinite coupling. We can then 
examine the NFL properties of the two channel Kondo 
model in a strong-coupling expansion. We can also ap- 
ply semi-classical methods to examine the model both in 
weak coupling and in a lattice generalization. 

Many aspects of the two channel Kondo model 
have been characterized by a combination of Bethe 
ansatz methods-*-^'-*-^, conformal field theory-*'*, numerical 
diagonalization^^'^^ and bosonization.^'*'^^ These meth- 
ods show that the NFL ground state is characterized by 

o Logarithmic temperature dependence of the specific 
heat Cv oc TlnT and magnetic susceptibility x 
InT. 

o A remarkable zero point entropy of magnitude | In 2. 

Emery and Kivelson's study^** of a solvable anisotropic 
limit of the two-channel Kondo model identifies this frac- 
tional zero point entropy with a unique fermionic zero- 
mode: a "Majorana" fermion. This localized mode de- 
velops below the Kondo temperature ultimately decou- 
pling from the conduction sea at the NFL fixed point. 
Sengupta and Georges^^ have shown that the residual 
low-energy coupling of this fermionic zero mode to the 
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conduction sea creates the logarithmic temperature de- 
pendence of the susceptibility and specific heat. Our 
work complements these discoveries and provides a sim- 
ple Hamiltonian formalism for studying them in the spin- 
isotropic limit. 

We begin, in Section II by compactifying the two chan- 
nel Kondo model, mapping the spin degrees of freedom 
of each channel onto the spin and isospin (charge) de- 
grees of freedom of a single-flavored electron fluid. In 
section III we make a strong-coupling expansion of the 
non-Fermi fixed point, applying methods previously ap- 
plied by Nozieres to the one channel Kondo model. In the 
strong-coupling Hamiltonian, non-Fermi liquid behavior 
is associated with the development of a sharp three-hody 
resonance in the conduction electron propagators. In Sec- 
tion IV we develop a versatile functional integral tech- 
nique for studying this model in the regime of weak cou- 
pling. We first use a mean field approximation to find the 
apparent order parameter, then we show that Z2 fluctua- 
tions of this order parameter restore the local symmetry 
and lead to the same non-Fermi liquid regime that was 
obtained in the strong coupling limit. Finally, in Section 
V we use the methods developed for an impurity model 
to formulate and solve its lattice generalization. Of par- 
ticular interest here, is the competition between "inco- 
herent" non-Fermi liquid behavior and coherent three- 
body bound-state motion. Section VI summarizes and 
concludes the results of these studies. 



II. MODEL 

Our starting model is the two-channel Kondo model, 
which we write in a one-dimensional form 

H' = it Y^[a{jn + l)ax^(n) - H.c] 



nXcr 



+ J[ai(0) + ^2(0)] • Sd 



(2) 



where A = 1,2 labels two independent one dimensional 
conducting chains and (J\{0) = aJ^(J(0)(Ta/3aA/3(0) is their 
spin density at the origin. In this model, prior to coupling 
to the local moment at the origin the spin degrees of 
freedom of the two chains are completely independent. 

Motivated by the observation that in the continuum 
limit, the spin and charge degrees of freedom of a sin- 
gle chain behave as two decoupled spin- 1/2 degrees of 
freedom, we now write down a compactified two-channel 
Kondo model where the local moment is coupled to the 
spin and isospin degrees of freedom of a single chain: 

H = it J2[ci{n + l)c,(n) - H.c] 

n 

+ J[a(0)+f(0)]-Sd. (3) 



Here 



is the spin density at position x and 

f(x) = cl(x)aai3Ci3(x) 

is the isospin density, whose diagonal and off-diagonal 
components respectively describe the charge and pair 
density at site x. The tilde on the electron operators 
denotes a Nambu spinor 



(4) 



We link the two models by making an identification 
of the low-energy spin modes in (2) with the low energy 
spin and isospin modes of (3). 



ai(x) 



a(x) 
f{x) 



(5) 



The dynamical equivalence of these two sets of opera- 
tors is established by comparing their long-distance, long- 
time correlators at J = 0. At J = 0, in a particle-hole 
symmetric band, the correlators of isospin and spin are 
identical. Furthermore, since the low energy spin and 
charge degrees of freedom are decoupled, the joint low en- 
ergy correlators of isospin t(0) and spin 't(O) completely 
factorize, thereby emulating the joint correlations of the 
spin density in the original two channel Kondo model. 
More precisely, 

(a(l)a(2)...a(n)f(l')f(2')...r>'))/=o 
= (ai(l)ai(2)...ai(n))(a2(l')^2(2')...'?(r'))j=o (6) 

if all time scales are longer than the band-width (i.e. all 
\Ti—Ti: I ^ 1/^)- (See appendix for a more detailed discus- 
sion). If we now turn on the interaction, this one-to-one 
correspondence between the correlation functions ensures 
that the models are equivalent in the continuum limit, 
to all orders in perturbation theory. By explicitly sep- 
arating the spin degrees of freedom in the compactified 
model, when we now use a lattice cut-off procedure we 
avoid developing an unrenormalizable coupling between 
spin and charge in the strong-coupling limit. 

It is convenient to rewrite the conduction electron 
spinor in terms of its real and imaginary components \f 
as follows 



c(n) = -^(*°(n) + i*(n).a) (^^. 



In terms of these components, (see Appendix B) 



a(n) + f(n) = -i\I>(n) x *(n). 



(7) 



(8) 



and 



a=3 



(4(n + l)c,(n) - H.c.) = ^ *<'(n + l)*<'(n). (9) 



a(x) = cl(x)aai3Ci3(x) 
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The Hamiltonian is seen to decouple into a "scalar" and 
"vector" part 



H = H., 



H,,. 



(10) 



where 



ff,e = it^*°(n + l)*°(n), (11) 

n 

H^=itJ2 *(" + 1) • - iJSd ■ (*(0) X *(0)) . (12) 



The "scalar" part (11) describes a freely propagating 
fermion mode whereas the "vector" part (12) describes 
a local moment interacting with a sea of vector Majo- 
rana fermions. This simple result shows that the low en- 
ergy spin excitations in the the two channels fuse to form 
a band of free singlet and interacting vector fermions. 
Schematically: 



spinon [1] (8) spinon [2] — > scalar ® vector. 



(13) 



From (8), we see that the joint spin-density of the two 
chains is carried by these vector excitations. 

The universal scaling properties of this model are iden- 
tical to those of the two channel Kondo model. The 
weak coupling scaling equation for the interaction con- 
stant g = pj (where p is the density of conduction elec- 
tron states at the Fermi energy) has a form 



(3(g) = 2g' - + 0[g''] 



(14) 



On general grounds, the first two terms in the weak 
coupling expansion of the beta function are universal. -"-^ 
Higher order terms which determine the location of the 
NFL fixed point can not be defined in a universal way, so 
the precise location of the NFL fixed point depends on the 
cut-off procedure. The original and compactified mod- 
els correspond to different regularization schemes with 
running coupling constants gi(^) and 32(0 that differ in 
strong-coupling due to the higher order terms in the beta 
functions. Their coupling constants coincide at weak cou- 
pling. 
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G[x] 



G[gi 



0\ 



(15) 
(16) 



but the higher order terms of G[x] depend on the cut-off 
procedure. In the compactified model there is no unsta- 
ble "over-screened" fixed point. Though the long-time 
correlation functions of 5 = |[(t(0) + t(0)] mimic the 

correlations of the operator S' = |['Ti(0) + 'T2(0)], in the 
two-channel problem 5 is a spin- 1/2 operator obeying 
usual identities: 
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(17) 
(18) 



On a lattice, the spin and isospin can not simultaneously 
be finite at one site, so it is not possible to produce 
a spin 1 composite at the origin to over-screen the lo- 
cal moment. Consequently, the unstable "over-screened" 
fixed point vanishes from the renormalization flows. In- 
stead, as we shall show, the non- Fermi liquid behavior ap- 
pears as a stable fixed point at infinite coupling g2 = oo. 
Qualitatively, if g\ is the zero of the beta function for 
the original two-channel Kondo model, by changing the 
cut-off procedure we have mapped this point to infinity 
9*2 = G[gl] = oo. 



III. STRONG COUPLING LIMIT 

In this section we examine the strong coupling limit of 
the compactified two-channel Kondo model. For simplic- 
ity we assume that the local moment is located at the 
origin of a chain: 



= it ^ *(n + 1) • *(n) - iJSd ■ ( *(0) X *(0) ) . (19) 



We imagine that we carry out a Wilson scaling 
procedure, keeping track of the corrections to the 
free energy and the evolution of the Hamiltonian as we 
integrate out the high energy degrees of freedom and 
move away from weak coupling. On general grounds, 
the qualitative physics described by the Hamiltonian at 
each point on the scaling trajectory will not change. If 
the strong-coupling fixed point is stable we can obtain a 
simple description of the essential physics of the model 
at weak coupling. 

We begin by examining the character of the strong 
coupling fixed point. At J = c» we may neglect all sites 
except n = 0: 



-iJSd • (* X *) 



(20) 



The eigenstates of the strong coupling Hamiltonian Hgtr 
correspond to singlet and triplet states formed between 
the local moment and either the spin or isospin of con- 
duction electron at the origin. Thus, there are two sin- 
glet and two triplet states, separated by energy 2J, as 
shown in Fig. 1. This two-fold duplication of energy lev- 
els results from a local supersymmetry: there are two real 
Fermi operators, \f and 

$ = -2i*(^)*(^)*(^), 

which commute with the Hamiltonian Hoo (20). The 
complex fermion Q = (\f — i$) /-\/2 transforms one state 
into its degenerate partner. Unlike bosonic degeneracies, 
transitions between the two singlet states is always ac- 
companied by emission/absorption of an odd number of 
fermions. 

To examine the effect of hybridization we first reinstate 
it between sites and 1, writing 
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S=1 



2J 

1' 



s=o 



n=1 



n=0,2 



Fig. 1.: Energy levels for the strong coupling 

limit. The singlet can form with a spin (n = 1) or 
"isospin" (n = 0,2). The two-fold degeneracy can be 
identified with the fermion ^ which commutes with the 
strong coupling Hamiltonian. 



H — Hoo + i, 



(21) 



where i = (1) • \f (0). The hopping term t causes tran- 
sitions amongst the singlet and triplet states at the ori- 
gin. Its matrix element between the singlet states is zero. 
We now use perturbation theory in t/J to eliminate the 
high energy states triplet states and form a low-energy 
effective Hamiltonian H*. This is done by making the 
canonical transformation: H — > H* = He~^ , where 



-2J 



t 



Pi - (H.c.) (22) 



is chosen to eliminate the off-diagonal terms between the 
singlet and triplet states. Here Pi and Pq respectively 
project into the subspaces with a triplet, or singlet con- 
figuration at the origin. The residual interaction induced 
in the singlet subspace is then 



Hint 



11 
"2J 



(2jy 



Pn. 



(23) 



Under the singlet projection, the first term in (23) van- 
ishes. The second term becomes 



Hi 



,$*(l)(l)*(2)(l)^(3)(l)^ 



(24) 



where a = 3t^/4J^. This term couples to the fermionic 
zero mode $ to a three-body composite of conduction 
electrons at site 1. In the more conventional one channel 
Kondo model, this term is absent and the leading term 
in the Hamiltonian is a benign Fermi liquid interaction 
proportional to f^/J^. The three-body resonance feature 
of the two-channel fixed point fundamentally transforms 
its physics. Our final fixed point Hamiltonian for the 
spin-degrees of freedom of the two-channel Kondo model 
is thus 



11=1 

+ a$*(i)(l)*(2)(l)*(3)(l). 



(25) 



where the site at the origin is explicitly excluded. 

Let us now consider the perturbative corrections that 
result from coupling the itinerant band to the zero 
mode. We shall show that these corrections are for- 
mally in the renormalization group sense, but that 
they produce strong scattering effects on the conduc- 
tion electrons. We construct the perturbation theory 
using Feynman diagrams containing the local propaga- 
tors of the itinerant conduction electrons and the three- 
body bound state. Both propagators describe Majorana 
fermions, and are represented by lines without arrows. 
We represent the conduction electron local propagator is 
{^-{l,T)^\l,Q))=Sai,G{T) as 



(1,0) 



(1,t) = SabGir), 



where 



PttT 



sin(T7rT) 



(26) 



(27) 



Here T is the temperature, eg = — 2tsinfc and p = 2^vf 
is the density of states for a Fermi velocity vp = 2t. 
The bare three-body bound-state propagator G^(t) = 
($(t)$(0))o is 







where 



G%{t) 



-TV 



:sgn(T), 



(28) 



(29) 



At T = 0, G(t) ~ so the scaling dimension of the con- 
duction electron operator \f(l) is d$ = |, whilst the 
fermionic zero mode is dimensionless: dq, = 0. The 
dimension of the interaction term in the action Si„t oc 
/ dT^*^*^*^ is consequently 



^int 



3d* = -1/2. 



The interaction is formally "irrelevant" , because it scales 
as [t"-*-/^] = [w'''''^]. From a practical standpoint, this 
means that each additional power of associated with 
vertex corrections to the skeleton diagrams introduces 
an additional power of frequency into the diagram. The 
vertex corrections to the bare interaction vertex between 
three conduction electrons and the zero mode, shown in 
Fig. 2 are thus proportional to a'^w and can be completely 
neglected at low energies. 

The polarizability of the electron fluid leads to two 
leading corrections to the interactions between the local- 
ized and itinerant fermions, as shown in Fig. 3. The ver- 
tex induced between the localized Majorana modes and 
the itinerant electrons (shown in Fig. 3a) is a'^Xeei'''), 
where Xee(T) is the local spin-susceptibility of the con- 
duction electrons. At long times, Xee(T) becomes small: 
XeeiT) OC which implies a frequency dependence 

Xeei^n) cx \i'„\ showing that this term is irrelevant at 
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Fig. 2.: Vertex corrections to the bare inter- 
action. Solid lines denote the electron propagator 
G(t). Dashed lines denote the three-body propaga- 
tor G^(t) = |sgn(T). The vertex corrections scale as 



5CeeW~0(1/T2)^ 101X33 (CO) ~N 
(b) 

1/T 




Xe^(T) ~ 0(1 /T) ^ImXe^ («) ~ In co 

Fig. 3.: Renormalization of interactions (a) Irrele- 
vant susceptibility XeeiT) oc (b) marginal susceptibil- 
ity between conduction electron and three-body bound 
state Xe^{T) oc K 



low frequencies. The second vertex shown in Fig. 3b de- 
pends on the joint-susceptibility between the conduction 
electrons and the three-body bound-state Xe#- From a 
scaling argument, Xe^(T) oc 1/|t| which implies an inter- 
action g^Xe^(^) oc In |a;|. This term generates a marginal 
interaction between the conduction electrons that is re- 
sponsible for departures from Fermi liquid behavior. 

Another way to consider the effect of these interac- 
tions is to evaluate the self energy parts that they gener- 
ate. Only two diagrams are possible in the second order 
(Fig. 4). The term shown in Fig. 4a renormalizes the 
propagator of the fermionic zero mode. 



G^(w) G$(w) = 



w — S$(a;) ' 



(30) 



Since S$ (t) oc I/t^, ImS$(w) oc w at low frequencies. 
This means that the sharp pole in the zero-mode propa- 
gator is preserved to all orders in the scattering processes 



-Im[G#(a;)] = Z6{w) + (background) 



(31) 



where = [1 — (5tjS$(a;)]|o. In other words, a renor- 
malized component of the three-body operator commutes 
with the fixed point Hamiltonian. By contrast S*(t) oc 
and hence /mS$(a;) oc |a;|. 



S$(a;) oc wlnw. 



(32) 



a signature of a marginal Fermi liquid. The stability of 
the fermionic zero mode guarantees that this feature is 
not removed by higher order processes. A consequence 



is that the itinerant spinous described by the conduction 
sea develop a sharp, zero energy pole in their three-hody 
propagator. 

The explicit form of the off-diagonal susceptibility Xe# 
in the frequency domain is of particular interest. We 
write 



Xe^(T)S'^' = ( [*«(!, t)$(t)] [*"(!, 0)$(0)]) 
so that 



Xe^(Vn) 



dTG{T)G%{T)e'--\ 



(33) 



(34) 



Writing this as a Matsubara sum 

sgn(w™) 



|w™|<Z» 

^ln^-*fl 



27rT 



'2 2'kT' 



(35) 



where \f (x) denotes the digamma function.^® Making the 
analytic continuation gives 



P 



2'kT ^2 2'kT' 
The imaginary part of this susceptibility is 
"1 



x'Uiv) = P 



„ tanhf — ) 
2 ^T^ 



(36) 



(37) 
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(a) Liz) ~ 0(1 /t3) =^ ImZ(co) ~ co' 




(b) E(x) ~ 0(1 /i^ ) =4 i:(co)~co In CO 





(c) r(coJ = -iAsgn(ccf, ) 

Fig. 4.: Self energy corrections (a) Irrelevant correc- 
tion to three-body self energy. Since ImT,{w) ~ w^, this 
term preserves the sharp zero-energy pole in the three- 
body propagator, (b) Marginal self energy of conduc- 
tion electrons induced by scattering off sharp three-body 
bound-state. /mS(a;) ~ |a;|. (c) "Hartree" self energy 
that gives Majorana bound-state a finite width A cx 
in a magnetic field B. 



T ' 

constant. 



(ly « T) 
» T) 



(38) 



The "marginal" form^ of this susceptibility is a direct 
consequence of scattering off the perfectly sharp three- 
body bound-state. 

Suppose we introduce a magnetic field by adding the 
term of a magnetic field that couples to both conduction 
electrons and the local moment, 



HB = B.[S,-ii/2)J2'^j 



(39) 



i=o 



The splitting between the singlet and triplet states at 
site is unaltered up to terms of order 0(B^ /J), so the 
residual coupling of the magnetic field in the fixed point 
Hamiltonian only involves the electrons at sites x > 



(40) 



These terms introduce an off-diagonal contribution to the 
conduction electron propagators 

(V-^l, t)V'6(1, 0)) = Go{T)8ai + ie^fteSepe-^l^lsgnr, (41) 

where p = 1I{2'kvf) is the density of states and D the 
bandwidth. The second term in the propagator gives rise 
to a static magnetization 



--(V(l)xV(l)) = Bp, 



(42) 



and the zero-mode self-energy now acquires a "Hartree" 
self energy, as shown in Fig. 4. This term breaks the 
fermionic degeneracy and broadens the sharp resonance 
into a Lorentzian of width A = ■jrB^(ap)^ p. The field- 
dependent energy scale will act as an effective Fermi tem- 
perature 



Tf ~ A(B) (X B^. 



(43) 



X^(v) 



5Ceci,(v) 



XJv) 




Fig. 5.: Impurity susceptibility. Showing how the 
magnetic field couples to the marginal susceptibility. 

One of the most interesting physical quantities to cal- 
culate is the dynamical magnetic susceptibility of the im- 
purity. This is given by the diagram shown in Fig. 5, 
which leads to 



(44) 



where Xci^^) is the uniform magnetic susceptibility of the 
conduction band. At energies \i>\ << t, Xci^^) = P, so 
that the impurity spin susceptibility is directly propor- 
tional to the marginal susceptibility Xe^i^^) 



Xspin(v) = (apfxe^(v)- 



(45) 



From this result, we may also deduce the the static sus- 
ceptibility 



Xspin = p(ap)Hn 



(46) 



where A = £)/(27re*( = )). 

We may calculate the dominant non-Fermi liquid cor- 
rections to the free energy within second order perturba- 
tion theory. These calculations essentially duplicate the 
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earlier work by Sengupta and Georges. In a weak mag- 
netic field B << T , there are two leading contributions 
to the free energy ( Fig. 6.), a zero field part (Fig. 6(a)) 

F,iT) = ^a' l\Gir)r[^sgnr]dr 
>(«p7rT)2] f 



7r/2 ^ 

dx 

J Jt/a 
./n ^ 



• 3 

sm X 



and a finite field part (Fig. 6(b)) 
(aM)2 



F2{T) 



piapBf 



L 



p{apBf 



In 



7r/2 ^ 

dx— 

x/A smx 

A" 



(47) 



(M = pB) 



(48) 



The bulk susceptibility and linear specific heat of the 
original two channel Kondo model are 



Xo 



2Np 



(51) 



where TV" is the number of sites in the lattice, giving a 
Wilson ratio 



W 



Xi/ Xo 

alio 



(52) 



At fields large compared with the temperature, the Ma- 
jorana zero mode behaves as a resonant level model with 
a finite width A cx B^. In this case, the repeated scat- 
tering off the resonance must be summed to all orders. 
This gives rise to a resonant-level contribution to the free 
energy 



F[A] 



2 ^ 



(53) 



(a) 



F,(T) : 




(b) 



F„(T) : 



Fig. 6.: Leading dangerous correction to free en- 
ergy, (a) Zero field component: leading correction from 
scattering off the three-body bound-state; (b) Finite field 
component. 

Adding these two contributions to the free energy of 
the asymptotically decoupled zero mode, Fq = Eo — '^lnT 
gives 



-InT ■ 



p(apy 



[{irTf + 2B^] In 
(B « T) 



(49) 



This leads to logarithmic contributions to the suscepti- 
bility and linear specific heat {j{T) = CviT)/T) 




(50) 



where A„ = sgn(u!„)0(D — |a;„|)A. Using the result 



In 



r(z + k) 



T(z) 



(54) 



0<n<k 

we can carry out the Matsubara sum to obtain 



F[A] = Tin 



V2^ 



^ln(^ 
27r V 2'kT 



(55) 



provided (A(B) << T). The function 

.r(| 



dF 
dT 



-In- 



27r 



A^ 
27r 



(56) 



describes the field-induced quenching of the |/n2 entropy 
associated with the three-body bound-state. 

To complete our discussion of the non-Fermi liquid 
fixed point, we should like to consider the entropy associ- 
ated with the fermionic mode. In the single site (J = c») 
limit there are two real fermions, \f and $ that com- 
mute with the Hamiltonian, which gives an entropy In 2. 
Once the coupling to the remainder of conduction sea is 
included, the scalar fermion \f becomes a part of the 
decoupled Fermi band, so half the fermionic degeneracy 
is removed. To further elucidate the nature of the "frac- 
tional" degeneracy it is instructive to consider two "two 
channel" Kondo impurities located at opposite ends of 
the conduction chain (Fig. 7). In such a situation there 
are two fermionic zero modes $(0) and $(7V") associated 
with each end of the chain. Using perturbation theory 
we may deduce that these two modes are coupled by 



= -iA$(0)$(Ar), 



(57) 
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where \ = a J G^{T)dT and 

G(t) = 1,t)*(i)(1,0)) 



2i 



N 



(vftY 



0{l/N^). 



(58) 



Thus A ~ /(vpN^). The interaction (57) leads to a 
splitting of two zero modes. To show this, it is convenient 
to form a complex fermion ^ = ■^[^(0) — i^{N)], then 



Hi 



(59) 



indicating that as TV" ^ c» the ground state of this sys- 
tem will develop a two fold degeneracy and, hence, a zero 
point entropy S = In 2. This entropy is completely de- 
localized, and for this reason we must associate a | In 2 
entropy with the residual real Majorana zero mode $(0). 
The reasoning employed here is strongly reminiscent of 
arguments used to assign a fractional charge to a soliton 
in a polyacetylene chain. 



<1.{0) 



*(N) 



N-2 N-1 




; 1/N^ 



S = ln2E (1/2)ln2 per impurity. 



To develop this generalization, we introduce a Majo- 
rana representation of the local moments^** 



--7] X 7], 



(60) 



where the operators f) are real and satisfy canonical anti- 
commutation rules {77", 7]^} = Sab- In terms of these op- 
erators the Hamiltonian can be written 



H 



ipin + 1) • il>(n) 



J(i?7 • *(0))2 



(61) 



We now generalize the number of components of rj" 
and Vf" from 3 to TV: a = 1 . . .N. Then the generalized 
Hamiltonian 



H 



ipin + 1) • il>(n) 



^(i,7•*(0))^ (62) 



has a global 0(N) invariance 



(63) 
(64) 



where Rab is an TV" dimensional orthogonal matrix. We 
have rescaled J —f J /2N so that with J fixed, the Hamil- 
tonian grows extensively with TV, establishing a well de- 
fined mean theory in the large TV limit. 

The interaction term can be conveniently rewritten in 
terms of 0(N) spins 



Sf = -i[v\v% 



(65) 
(66) 



whereupon 



Fig. 7.: Fractional entropy. A two channel Kondo 
chain with TV + 1 sites. For finite chain length the two 
Majorana zero-modes are coupled, giving rise to a split- 
ting of the ground state of order 0(1/TV^). In the limit of 
TV ^ c», every state in the Hilbert space is two-fold de- 
generate, leading to a "delocalized entropy" of |ln2 per 
local moment. 



IV. WEAK COUPLING LIMIT 

In this Section we solve a generalization of the spin- 
charge Kondo model (12) in which the number of com- 
ponents in the vector Majorana is generalized from three 
to TV. We then formulate the model as a path integral 
and show that the fluctuations about the large TV limit 
(i.e. about the saddle point solution) recover the key as- 
pect of the strong coupling limit, i.e. they reproduce the 
effective low energy Hamiltonian (25). In contrast to the 
previous section, this method can be applied to the weak 
coupling properties of the model. 



Hint 



2TV 



SfSf(O). 



(67) 



The 0(N) Kondo model possesses a local Z2 invariance 
under transformation 



V 



-V- 



(68) 



The explicit form of the Z2 operator depends on whether 
TV is even or odd. For the case of TV even 



■V 



(69) 



The case of odd TV requires special consideration. To 
construct the Hilbert space which represents the Majo- 
rana fermions rj^jT]^ . . .r]^ , we may pair TV — 1 fermion 
rf' . . . 7]^ into complex fermions 



ci 



V2 ' 



2, 4. ..TV- 1. 



(70) 



To deal with the remaining rj^ we are obliged to introduce 
an additional "ghost" fermion field $ such that if rj^ = 



co+c 



■yj^, $ = — . With this representation, for TV odd: 
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iV+1 1 



N 



(71) 



Since $ commutes with the Hamiltonian, this also implies 
that for N odd, the fermion operator 



N 



(72) 



commutes with Hamiltonian. This can be also verified 
directly. Thus, for odd TV", the Z2 invariance manifests 
itself as a supersymmetry 



[2, H] = Q. 



(73) 



We now formulate the partition function as a functional 
integral 



S= / -(775,77 + V*(n) •5,*(n)) +ff 
Jo ^ 



(74) 



dr. (75) 



To proceed, we make a Hubbard Stratonovich transfor- 
mation on the interaction term 



J 

"2N 



irf ■ \f 



iy(T)(77- *) 



2J 



(76) 



where the V(t) is a fluctuating real field. In the large TV" 
limit, fluctuations of this field are small, and at TV = c» 
the thermodynamics is determined by the mean-field sad- 
dle point V = jf{irj ■ ^(0)). This saddle point solution 
breaks the local Z2 invariance, and at finite TV fluctua- 
tions in the sign of V{t) play an important role in restor- 
ing the local symmetry. 

The mean field Hamiltonian 



Tvy2 



^rMF(I^) = ^efc*-fc*fc+I^[i»?-*(0)] , 
4k = TV-i/2^*(j)e-'^^^ 



(77) 



describes an TV-fold degenerate band hybridizing with a 
resonant level. At this level of approximation, the prop- 
erties of the Fermi system are those of non-interacting 
resonant level of width A = irpV^. Using the resonant 
level free energy obtained in (55), we find 



Fmf = -T Tre"''^"^ = N{F[A] ■ 



A 



27rJp 



Minimizing this with respect to A, we obtain 



V'(; 



27r 



) - In 



27r 



where 



D exp 



(78) 



(79) 



(80) 



is the "Kondo temperature" which sets the characteristic 
energy scale of the mean-field properties. At low temper- 
atures, A — > Tk- 

Let us now discuss the fluctuations about the mean 
field theory. In the functional integral 



/ 



X>*X>77X>ye" 



where 



I {vdrV + *(«)5r*(n)) + Hmf(V) 



(81) 



dr, (82) 



we may identify two classes of fluctuations: (i) small 
Gaussian fluctuations around the saddle point solution 
and (ii) rare tunneling events when the sign of the real 
field V reverses ("kinks") as illustrated in Fig. 8. 

In the large TV limit, the action associated with an 
order parameter "kink" scales as TV, and hence the fre- 
quency of these fluctuations scales as r(TV) cx exp(— TV). 
The small amplitude Gaussian fluctuations of frequency 
w ^ r will be unaffected by the rare large amplitude 
kinks, thus in the large TV limit we may consider the 
contributions to the functional integral from these fluc- 
tuations independently. 

Consider first the effect of the Gaussian fluctuations. 
Expanding the action around the mean field theory we 
find 



S = S 



TV 



MFT 



J2 n(z^„).5y(iz.„),5y(-z.„) + oisv^), (83) 



where H(j^„) is given by the fermion bubble 



IL(iiy„) 



TV 



((i77-*)_^„(i77-*)^J 



(84) 



Since Tk is the only energy scale in the problem, at low 
temperatures H(j^„) has the functional form 



(85) 



Though we can calculate f(x), its detailed form is not 
of great importance. Since there is no continuous broken 
symmetry associated with the mean-field theory, these 
Gaussian fluctuations are gap-full at all frequencies and 
f(x) is a positive analytic function of x. Gaussian fluc- 
tuations therefore generate small 1/TV corrections to the 
dynamics that are analytic functions of the frequency, so 
they do not give rise to non-Fermi liquid behavior. 

Next we consider the effect of the large discrete fluctu- 
ations in the sign of V . Let us divide the trajectory V(t) 
into periods of time when V = ±Vo is almost constant, 
separated by brief periods when V(t) is changing sign. 
The time scale of the sign change can be estimated by 
the rate at which V(t) approaches the stationary value 
±Vo after the transition has occurred. This rate can be 
estimated from (84): I/tq ~ Tk- 
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Free Energy 




Fig. 8.: Kink configuration. Configuration with four 
kinks in the order parameter y at a given site. 

We should like to compute the action associated with 
2n widely spaced sign changes in V(t). We anticipate 
that there will be fixed amount of action associated with 
each transition, but there will be also Berry phase con- 
tributions that lead to long time interactions between 
kinks. 

To calculate the overall Berry phase contributions to 
a tunneling trajectory we may use the "sudden approxi- 
mation" in which each tunneling is considered to be in- 
stantaneous. The contribution to the path integral from 
a given tunneling trajectory of 2n kinks is then 

T^"Tr [^e-^] , (86) 

where F is the amplitude of a single tunneling event and 

S= H{Vo)dT- I H{-Vo)dT...- ( H{Vo)dT. 

(87) 

If we take advantage of the identity 



Te 



■r 

•J T' 



H(-Vo)dT 



PTe 



H{Vo)dT . 



(88) 



we may write the contribution to the path integral from 
a given tunneling trajectory in the form 



-Tr 



Tnp(r,)a.e-r^(^°)<'^ 



i=i 



, (89) 



where the trace includes the trace over a pseudo-spin 
variable a = sgn[y] that keeps track of the sign of Vo, 
and the normalization factor of | takes care of the over- 
completeness introduced by introducing both the Majo- 
rana spin representation and the additional ghost field 



The operator converts pseudo-spin up to pseudo-spin 
down and vice versa, so Am [{Ti}] vanishes for m odd, and 
includes the contributions from the two trajectories with 
flips at times {tj}, starting out with either 1^(0) = Vo 
or V(0) = —Vo- The total contribution to the functional 
integral from all such trajectories is then given by 



A„ = — ( dTi. . .dT„A„[{Ti}]. 
n\ J 



(90) 



We may now recognize the sum of all A„ as the per- 
turbation expansion for the following time-ordered expo- 
nential 



-Tr 



(91) 



To estimate F, we may set it equal to a product of the 
attempt rate (~ Tk) and the success rate of the tunneling 
process, the latter which becomes exponentially small for 
large TV". To find the success rate with exponential accu- 
racy we need to find the optimal trajectory leading from 
one minima to another which minimizes the action. We 
shall not do the explicit calculations here, instead we es- 
timate it as an action associated with a very fast change 
of the sign. The exponential of this action is the overlap 
between the wave functions of the fermions in the ground 
states before and after the sign change, so 



(92) 



We estimate (+|-) ~ exp(-7V"/2). 

We may rotate in the pseudo-spin space to transform 
o'x —^o'z- Since cr^ commutes with the Hamiltonian, we 
may choose one sign of the pseudospin. The factor of 2 
in the partition function that this gives rise to, cancels 
the normalization in (91). The effective Hamiltonian for 
low energy fluctuations then becomes 



Hn = H(Vo) + TP. 



For even TV", 



HN = H{Vo) + fv^...V 



N 



(93) 



(94) 



describes a Fermi liquid, with a weak interaction amongst 
the hybridized Fermi fields rj. However, for odd N, the 
effective interaction couples the zero energy Majorana 
fermion which is un-hybridized, to a composite of TV 
hybridized Majoranas r]^ . . . ■q^ : 



(95) 



This result is qualitatively similar to that obtained in the 
strong coupling analysis of Section III, but it is now seen 
to be a feature of all 0(N) Kondo models with odd TV. 
By repeating the perturbation arguments of the previous 
section, we now see that all odd TV, there is a decoupled 
Majorana zero mode with a | In 2 zero point entropy. 
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Although the effective Hamiltonian (95) has the same 
qualitative form for all odd N as the effective Hamilto- 
nian (25) of the strong coupling limit, the large effect 
of the Majorana zero mode on conduction electrons is 
specific to TV = 3. Consider, for example, the electron 
self-energy produced by weak scattering off the induced 
on-site interaction. The corresponding diagram is shown 
in Fig. 9. Since, the self energy scales as 1/t^~^, this 
implies 



/mS(a;) ~ w 



N-2 



(96) 



Marginal behavior only occurs for TV = 3. The correction 
to the free energy for odd N 



InT, 



TV = 3 
TV > 3. 



(97) 



Thus non-Fermi liquid behavior only develops for the 
physical 0(3) Kondo model and the unusual fermionic 
zero mode produces only weak corrections to Fermi liq- 
uid behavior at TV > 3. 




(N-l) 



thus be to cut-off the non- Fermi liquid behavior by elim- 
inating spin-charge decoupling. 

In this section we consider the effect of the lattice on 
the zero-mode behavior of the 0(3) Kondo model. Unlike 
the two-channel Kondo lattice, here the complications of 
these spin-charge coupling effects do not arise. We begin 
by writing down the "0(3)" Kondo lattice model 



H 



C jCj+a, + H.C 



'Ji ' 



(98) 



where j represents the site on a three dimensional cubic 
lattice, 



(99) 



denote the spin and isospin density at site j and a; (I = 
1, 2, 3 denotes the three unit lattice vectors. By using the 
Majorana decoupling of previous sections: 







(100) 



where Q = {t: 12,1: 12,1: /i), this model can be re-written 



where the decoupled scalar degrees of freedom, described 

by 



(102) 



Fig. 9.: Conduction electron self-energy in 0(TV) 
Kondo impurity. S*(t) scales as jv-2 ? which implies 
ImS$(a;) cx w^"^ is marginally relevant when TV = 3. 



V. THE 0(N) LATTICE. 

In this section we apply the approach developed in Sec- 
tions III and IV to a lattice generalization of the two- 
channel Kondo model (2). The two complimentary meth- 
ods of the last two sections have demonstrated how the 
driving force for non-Fermi liquid behavior of the two- 
channel impurity model is the formation of a fermionic 
zero mode. How specific is this behavior to impurity 
models, and do we expect similar type of behavior to 
persist in lattice generalizations of these models ? 

For more than one two-channel Kondo impurity, back- 
scattering between spin-sites will couple "right" and 
"left" moving currents, so the kinematic spin-charge de- 
coupling present in the original model is lost. In a two 
channel Kondo lattice, the main effect of the lattice will 



may be neglected. By writing the local moments in a 
Majorana representation, Sj = —-^ffj x ffj, this model 
may be cast in a manifestly 0(3) symmetric form 



H 



Actually, both forms of this model are of interest. The 
model written in form (98) involves a coupling between 
pair, and spin degrees of freedom. Such a model ap- 
pears very naturally as the low-energy model for an odd- 
frequency^® superconductor with composite off-diagonal 
order. In the more symmetric form, (103), the model 
can be compactly generalized to an 0(TV) model, by ex- 
tending the number of components of the vector fermions 
from 3 to TV, as was done in the last section: 



J 

2TV 



By examining this model around the large TV limit, fol- 
lowing the methods of the last section, we can learn how 
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the fermionic zero-mode and the non- Fermi liquid behav- 
ior are modified in a lattice environment. 

Following the procedure introduced in Section IV 
we decouple the exchange interaction by the Hubbard- 
Stratonovich field V(x, t) to obtain the following action: 



J 



In the mean-field approximation the y-field acquires a 
non-zero average Vo = {V) . To find we need to solve the 
equation 



SF(Vo) _ NVo 



J 



+ i((77j.Vi)) = 0, 



(105) 



where the fermionic Green's function is calculated against 
the background of a fixed uniform V. The mean-field 
Green's function on this background 



G(w„,k) = 
is given by 

G(w„,k) 



{{rr)) {{rv")) 



-ej iV 
-iV iu>„ 



(106) 



(107) 



where eg = 2t^j_^ sinp;. The quasiparticle 
spectrum, determined from the eigenvalue equation 
det[G(£;(k),k)-i] = 0, is 



%A = Y + We'(k)/4 + X^o^ 



(A = ±l) (108) 



contains a small gap: 



m)\ 



(109) 



The mean-field free energy per unit cell is then 
NT 



[Vo] 



^ln[2cosh(/3%^/2)] 



k,A 



NVl 
2J ' 

Minimizing with respect to Vo yields 



yo = J^tanh(/3^j /2) 



(110) 



(111) 



Replacing the momentum integral by an energy integral 
near the Fermi energy — > P J de and performing the 
integral at T = yields 



Vo = 2pJVo In 



D 

\Vo\ 



(112) 



where D ^ 6t is the upper cut-off and p is the density of 
states at the Fermi level of itinerant electrons. Solving 
this equation for Vo with logarithmic accuracy we find 



Vo K, De 



(113) 



where Tk = D~'p^ , so that the gap Ag = Tk- 

Clearly, the Z2 symmetry implies that the mean-field 
free energy is independent of the sign of Vj at each site, 
and beyond mean-field theory, V^(t) will fluctuate be- 
tween Vo and —Vo- We can apply the technique devel- 
oped in Section IV to take into account the effect of these 
sign fluctuations. In particular, it was shown that this 
effect can be emulated by the integration over an addi- 
tional Majorana ghost field $: 



cx j Z)$exp - j dT\^<^dr 



(114) 



where the coupling constant g is physically identified with 
the inverse tunneling rate for a local change of sign in 
the y(n, t) field. As local processes are controlled by the 
energy scale Tk, we may estimate g ~ Tk- The scale 
Tk is exponential in the bare coupling constant J, so the 
result (114) is essentially non-perturbative. 

The fermionic variable $ represents low-lying excita- 
tions of the theory. In order to derive the effective ac- 
tion for these excitations, we need to integrate out high- 
energy degrees of freedom described by the gapfull modes 
rji, ipi (i = 1, 2, 3). The mode $ interacts with these gap- 
full modes via 



Hint 



j 



(115) 



The simplest diagram for the self-energy part of $ is 
shown in Fig. (10). It results in the dispersion 



S(fc)~(^sinfcO^. 



1=1 



(116) 



Thus the low energy properties of the two-channel Kondo 
lattice are determined by a single band of real fermions 
with a bandwidth of order of the Kondo temperature. 
Since the fermion is real, the total amount of entropy ac- 
cumulated in this band is In \/2 per site, coinciding with 
the residual entropy of the single impurity problem. The 
effective inter-site interaction induced by virtual three- 
body fluctuations removes the degeneracy of the single 
impurity problem and results in the coherent dispersion 
of the $-mode. 

At low temperatures the propagating $-mode now 
leads to a linear specific heat. Since the matrix elements 



12 




under the mapping 



Fig. 10.: Inter-site coupling of three body bound- 
state. The three-body bound-state can disperse from 
site to site via virtual fluctuations into the conduction 



of spin operators and $-states are zero, there are no low- 
energy spin-excitations. Remarkably, the two-channel 
Kondo lattice represents a spin paramagnet with a finite 
spectral gap 2Tk- Above the energy gap a single spin 
kink is a composite excitation consisting of two gapfull 
Majorana fermions. In this respect this model is similar 
to the spin- 1/2 Heisenberg chain. 

In the weak coupling limit of this model, the three- 
body bound-states develop coherence at a temperature 
9^ /Tk ~ Tk- Clearly then, the simplest lattice general- 
ization of the 0(3) model will not display any non- Fermi 
liquid like behavior. In the strong-coupling limit how- 
ever, there is a possibility of a large separation between 
the the scale /Tk and Tk, leading to a regime where 
the three-body bound-states are localized and conduction 
electron self-energies are marginal. This is an interesting 
issue that we return to in the final section. 



VI. CONCLUSION 

A goal of this paper was to explore the non- Fermi liq- 
uid properties of the two-channel Kondo model by using a 
"stripped-down" version of the this model that describes 
the decoupled spin degrees of freedom. The compactified 
model that arose from these considerations provides a 
new perspective on the two-channel Kondo model and 
opens up a new family of 0(N) Kondo models which 
generalize the essential physics beyond the realm of the 
original model. In this final section we discuss further 
applications of our methods and touch upon some of the 
interesting questions that arise in connection with the 
0(N) Kondo lattice model. 

The key observation of this paper was the equivalence 
of the two-channel impurity Kondo model 

H' = itJ2[a{jn+l)ax^(n)--E.c.] 

n 

+ J[ai(0) + ^2(0)] • 5,,, (117) 
with its compactified counter-part 

H= i t^[4(n+l)c,(n)-H.c.] 

n 

+ J[a(0) + f(0)]-Sd, (118) 



ai(x) 



f(x), 



(119) 



that resulted from spin-charge decoupling. The universal 
spin physics of both models is the same, but by preserving 
spin-charge decoupling in the compactified model, the 
geometry of the scaling flows away from weak coupling 
was deformed, moving the location of the NFL fixed point 
from intermediate coupling in the old coupling constant, 
to infinite-coupling of the redefined coupling constant. 

Several other applications of this procedure have yet 
to be explored. One clear application is the use of 
the the new formulation in an actual numerical Wilson 
renormalization calculation. Wilson's original numeri- 
cal renormalization group approach to the one-channel 
Kondo model took great advantage of a strong-coupling 
expansion to connect the numerical results to the strong- 
coupling fixed point. This can not be done in the stan- 
dard formulation of the two-channel Kondo model, but 
it can now be done with the new formulation. This then 
opens the possibility for a direct numerical calculation 
of the marginal Fermi liquid spin-susceptibility of this 
model. 

Another simple application of this reformulation is to 
the channel-anisotropic two channel Kondo model, where 
the original interaction 



Hint = [Jl'?l(0) + J2'?2(0)] • Sd, 

may be compactified into 



Hi 



[Ji'?(0) + J2t(0)] • Sd. 



(120) 



(121) 



Bethe-ansatz work on this model suggests that the char- 
acter of the paramagnetic fixed point changes subtly 
in response to channel anisotropy.^-*- Unlike the origi- 
nal model, the strong-coupling limit of the compactified 
channel-anisotropic model is a stable fixed point. The 
"Fermi liquid" of spins that results from such a descrip- 
tion appears to very naturally described as a band of 
vector and scalar Majorana fermions. 

Perhaps the most fascinating aspect of the compacti- 
fied model concerns the appearance of three-body bound- 
states and their close link with marginal Fermi liquid be- 
havior. The original marginal Fermi liquid hypothesis 
was advanced as a phenomenological framework to de- 
scribe the normal state of cuprate superconductors.^ An 
essential feature of the marginal Fermi liquid ansatz is 
the appearance of a two particle susceptibility 



T ' 

constant, 



(w « T) 
(w » T) 



and a related electron self-energy 

S(w) ~ wln(max[w,T]). 



(122) 



(123) 
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It proves remarkably difficult to construct even a toy 
model that furnishes a marginal self-energy down to ar- 
bitrarily low energies. A scale-invariant self-energy of 
this form requires a massless excitation. Almost in- 
evitably, beyond leading orders in perturbation theory, 
a nominally massless excitation develops a mass. In an 
attempt to try to overcome this difficulty, Varma and 
Ruckenstein^^ tentatively suggested that the massless ex- 
citation responsible for marginal behavior might have its 
origins in a three-fermion bound state. Quite remarkably, 
the marginal Fermi liquid behavior associated with the 
0(3) Kondo model is directly related to the formation of 
a three-body bound state. The fixed point Hamiltonian 



a. 

rp^ 



(127) 



E* =it^ *(n+ 1) •*(n) 
11=1 

+ a$[*(i)(l)*(2)(i)^(3)(i)]^ 



(124) 



owes its marginal Fermi liquid properties to the coupling 
of a Majorana three-body bound-state to the continuum. 
The real character of the three-body operator $ = 
means that its square is a c-number = |{$,$'}, so 
that any "mass term" of the form m^$^ = m^/2 can 
be absorbed into the Ground-state energy, which pre- 
vents the development of a mass term associated with 
this excitation. A Majorana three-body bound-state in 
an impurity model is thus massless in to all orders in 
perturbation theory. This essential Majorana character 
of the three-body bound state thus ensures that: 

• The response function of the three-body bound-state 

retains a zero-energy pole, giving rise completely 
scale-invariant response functions. 

• The joint susceptibility Xe#(i^) of the electrons and 

the three-body bound-state is marginal, with no 
characteristic energy scale 



Xe#(w)" = p 



■ tanh 



T 

constant 



(w « T) 
(w » T) 



(125) 



The generalization of the strong-coupling Hamiltonian 
(124) to a lattice environment is particularly interesting: 

+ «^$(i)vl>(i)(i)vl>(2)(j)vi>(3)(j-). (126) 



providing a rather interesting toy model for detailed cal- 
culations of marginal Fermi liquid behavior in a lattice 
environment. 

Another suggestive aspect of this work this work con- 
cerns its possible links with odd-frequency pairing. ^^'^^ 
The semi-classical methods that become exact in the 
large- TV limit of the 0(N) Kondo model are almost iden- 
tical to those used to develop a mean-field theory for 
odd-frequency pairing in the Kondo lattice. Consider 
the family of Kondo lattice models 



where ^(j) and f(j) denote the spin and isospin den- 
sity at site j respectively. This type of model was re- 
cently suggested as a low-energy Hamiltonian for an odd- 
frequency superconductor,^^ where the A dependent term 
develops in response to the growth of composite order 
{o'a(j)Tb(j)) = Rabij) between spin and charge degrees 
of freedom. The lambda term can be considered anal- 
ogous to a Weiss field for odd-frequency pairing. The 
model described by .ff [1] is the 0(3) Kondo lattice model 
considered in this article; .ff [0] is the conventional Kondo 
lattice model. Remarkably, if we use a Majorana repre- 
sentation of the local moments, Sa = ^ Vj ^^id the 
conduction electrons, then the same mean-field ansatz 



(129) 



can be used to extend the semi-classical methods used 
in this article to all values of A. The mean-field theory 
contains a spin and charge gap for each value of A, and 
as A evolves, the state evolves continuously. For A = 1, 
this semi-classical theory corresponds to that described 
in this article. In the limit A ^ 0, the mean-field theory 
provided here corresponds to a broken-symmetry state 
with composite order and odd-frequency pairing. These 
results suggest that there may be a whole range of qual- 
itatively similar physics that evolves continuously from 
the 0(3) Kondo lattice model to the conventional Kondo 
lattice model. 
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This is a model with "pre-formed" three-body bound- 
states. Clearly, marginal Fermi liquid behavior will per- 
sist in this model down to temperatures T ~ T*, where 
the three-body bound-state begins to propagate coher- 
ently. This temperature T* is set by the size of the cou- 
pling constant 



APPENDIX A: Validity of the compactification 
procedure. 



Here, we discuss in more detail, the compactified for- 
mulation of the two-channel Kondo model. We begin 
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by linearizing the conduction bands around the Fermi 
energy to write the two channel Kondo model in the rel- 
ativistic form 



H 



M r . 

A = l L 



is the spin density at x, 

f(x) = cl(x)aai3Ci3(x) 



(130). 



This model describes M chains of right-moving electrons 
with group velocity A, each interacting with a single 
spin Sd at the origin. Each right moving electron spinor 
ipxix) represents an electron scattering off a Kondo im- 
purity in an s-wave scattering channel. The properties of 
this model are determined by the correlation functions 
of the total conduction electron spin density in the non- 
interacting system. For instance, consider the perturba- 
tion expansion of the partition function 



'7 ^ Til 

n = l 



(131) 



Here 



B„= / n'^^iW(n)...4''(^n))o 

•'° i=l 

X (T5«nri)...5^(r„))o, (132) 

where the subscript (0) denotes correlation functions in 
the non-interaction system and 



'^tW= E V'*a(^)<t<'Va(^), 



(133) 



A = 1,M 



is the total conduction electron spin density. In a one- 
dimensional problem, the spin densities behave as inde- 
pendent degrees of freedom; for the M channel prob- 
lem, the the Kinetic energy decouples into a sum of 



E 



A = 1,M 



where 



independent terms Hj-in 
Hx = vp : ^xil) • ^\(—<l) '■ is a bilinear in the Fourier 
transform of the spin density <Tx{q) = Efc V"' Afc'?V'Afc+g at 
momentum q, where ipxk denotes an electron of momen- 
tum k. The dynamics of these modes are set exclusively 
by their commutation algebras ("current algebra"). The 
spin density satisfies the current algebra 



2ieabccrl{q - q') + q6''^6q,q< ■ 



5AA' 



(134) 



For the case where M = 2, the spin dynamics of two 
spin channels may be mapped onto the spin and isospin 
dynamics of a single chain. Consider the single chain 
model where both spin and isospin interact with a local 
moment 



is the isospin density at x, written in terms of the Nambu 
spinor 



i>{x) 



il>^{x) 
i,\{x) 



(136) 



For one-dimensional electrons moving with linear disper- 
sion, spin and isospin form two independent current al- 
gebras 



K(g), a\-q')] = 2ieabc<T'{q - q') + q^''^'^. 



[r'-iq] 



A- 



2ieabcT'{q-q') + qS'"'S. 



K(g),rV)] = 0, 



(137) 



(138) 



(139) 



where the label q denotes the Fourier transform of the 
density at wave vector q. Since the electron Kinetic en- 
ergy can be expressed as the sum of two independent 
bilinears of charge and spin current densities, 

-i../(V--,vv.,).» 

= ..f 2 : a(q) ■ ^(-.i) + f(,j) ■ f(-<i) :, (140) 

9 

it follows that spin and isospin behave as two completely 
independent spin degrees of freedom with precisely the 
same dynamics as the spin densities of two independent 
chains. When interactions are introduced that couple 
to to (Ji(x) and a2(x), their effect will be identical in a 
compactified model where we have made the replacement 



ai(x) 



a(x), 
f(x). 



(141) 



For example, in the case of channel symmetry in the com- 
pactified model, the local moment couples to the total 
spin density T(x) = ^(x) + f(x). This operator obeys 
same current algebra as ^t(x) = ^i(x) + ^2(x) in the 
original two channel model, so both operators have pre- 
cisely the same correlation functions in the absence of 
interactions. Thus in the expansion of the partition func- 
tion of the compactified two channel model 



Z_ 



°" Til 

^ n! 

11=1 



(142) 



where 



H = —ivp 
where 



dxi^K'^i'sdx + JSd-[<?{0) + f{0)], (135) C„= T[dTj{TT^'{Ti)...T^'^{T„))o 



X (T5-(ri)...5-(r„))o, 



(143) 
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the correspondence of correlation functions guarantees 
that the expressions in (132) and (143) are equal: C„ = 
B„ , establishing the equivalence of the two partition func- 
tions. Similar perturbative arguments can be made for 
any response function involving spin degrees of freedom, 
even in the case of channel anisotropy. 

Of course, if we now introduce a cut-off into these mod- 
els by rewriting them as one-dimensional lattice models 
with band- width 2t, the equivalence of the two channel 
model with its compactified counterpart is only guaran- 
teed in the limit J /t « 1, at temperatures and fre- 
quencies « t. If we apply a Wilson renormalization 
procedure to the two models, then although the two cal- 
culations will predict the same low temperature proper- 
ties, the details of renormalization flows will differ be- 
yond weak coupling. General scaling principles ensure 
that the the universal low energy dynamics predicted by 
the Wilson procedure will match up for the two models. 



APPENDIX B: Simple manipulations of Majo- 
rana Fermions 



Here we present some of the technical manipulations 
required to express the spin and isospin density at the 
origin in terms of the Majorana fermion operators. We 
begin with the expression for the Fermi field at the origin 







or more explicitly 



1 

7! 



^1 _ 

_ ^^0 



(144) 



(145) 



where (\f°,\f) are four real (or Majorana) components 
which satisfy 



(146) 



To deal simultaneously with the spin and isospin oper- 
ators, it is convenient to rewrite the Fermi spinor as a 
Balian-Werthamer four spinor, as follows 



( 



where 



/ » 

—i 
i 

V 



(147) 



(148) 



In terms of this four-spinor, the spin and isospin opera- 
tors are 



a = -C^aC, 



where 



(T = (T (8) 1 



(149) 



(150) 



and 

(Zi, Z2, Zs) 



(151) 



When sandwiched between the spinors Z and Z'^ , these 
matrices give the following expectation values: 

Z^SZ = Z^fZ = 0, 
Z^cr^cr^Z = 26"^ 
Z^a^T^Z = -26"^ 
Z^cr^cr^cr^Z = 2ie''^\ 



Z^cr''cr''T''Z 



-2ie'' 



(152) 



With these results, we can expand the expressions in 
(149) to obtain 



-Z^ 



(*o - is. ■ *)o:''(*o + • *) 



(153) 



and 



(*o - iff • *)t'' (*o + iff • *) 



Combining these two results, we find 
(T + f = -i* X 



(154) 



(155) 
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